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CN , Abstract 

Ph , When two or more subsystems of a quantum system interact with each 

, ^ ' other they can become entangled. In this case the individual subsystems 

can no longer be described as pure quantum states. For systems with 
only 2 subsystems this entanglement can be described using the Schmidt 
decomposition. This selects a preferred orthonormal basis for express- 
ing the wavefunction and gives a measure of the degree of entanglement 
. present in the system. The extension of this to the more general case 

' of n subsystems is not yet known. We present a review of this process 

1^-^ , using the standard representation and apply this method to the geometric 

• algebra representation. This latter form has the advantage of suggesting 
\^ ' a generalisation to n subsystems. 

o ■ 

1 Introduction 

^ ^ Quantum entanglement in 2-particle systems is currently well understood (for 

a useful review, including an extensive list of references, see §). But the quan- 
' tum behaviour of many-particle systems is more complicated and less well un- 

^ , derstood, and it is these systems that are of interest experimentally. The main 

limitation to the theoretical understanding of such systems is that the techniques 
which have been developed to analyse 2-particle systems do not easily gener- 
alise. Geometric algebra has the advantage that the number of particles being 
^> , analysed dictates the size of the space but otherwise does not alter the analy- 

• sis used. In this way, results developed in simple cases (such as the 2-particle 
system) can be more easily generalised to the n-particle case. 

In this paper we focus on 2-state quantum systems in the cases of one and two 
particles. We start by reviewing the standard matrix-based approach to single 
and two-particle pure states. We describe the Schmidt Decomposition, which 
provides a measure of the degree of entanglement present in a given system. 
We then introduce the density matrix to describe both pure and mixed states 
in a unified manner. We next turn to an analysis of the same systems using 
the multiparticle spacetime algebra framework developed by Doran, Lasenby 
and Gull |^, |^ . As a simple application we review the properties of the spin 
singlet state, frequently encountered in discussions of the Bell inequalities and 
EPR-type experiments 
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2 Single-Particle Pure States 



If there is only one particle present in the system then the spinor, \ip) eTL, can 
always be written in the form 

|V')-co|0)+ci|l) (1) 

where |0) and |1) are some pair of orthogonal basis states and co and ci are 
complex coefficients. Alternatively, all information about the state of the par- 
ticle can be expressed in a polarisation (or spin) vector, P, whose components 
are given by 

^ (a,) = (2) 
where the are the Pauli matrices 

^1^(1 J), ^2-(° 7), a3=(^j \y (3) 

It follows from this definition that |P| = 1. In this way the spin state of 
the particle can be expressed graphically as a point on the 2-sphere, and any 
evolution of the state of the particle can be thought of as a rotation of the 
polarisation vector. In many applications this sphere is known as the Block 
sphere. 



3 2-Particle Systems 

Suppose that two particles are described by states belonging to individual 
Hilbert spaces Tii and 7^2. The joint Hilbert space for the interacting system 
is Til ® , consisting of complex superpositions of tensor products of states in 
the individual spaces. A basis for Tii (g) 7^2 is constructed by taking the tensor 
products of the basis vectors for Hi and 7^2. Therefore any pure state of the 
composite system, ji/i) £ Hi (^0.2, can be expressed as 

^ \j) = H c,,j|z,j). (4) 

i,j=OA ij=0,l 

If more than one of the Cij^s are non-zero then each subsystem is no longer in 
a pure state and the system is entangled. In order to quantify the degree of 
entanglement we re-express equation in the form 

\ip) = cos(a/2)|0', 0') -I- sin(a/2)|l', l'), (5) 

which can always be done via a suitable change of basis. Here, |0') and |1') 
are orthonormal vectors in Hi and H2 and < a < tt/2 (so that cos(a/2) > 
sin(a/2)). Then |0',0') can be thought of as the separable state 'closest' to lip) 
and a as the degree of entanglement present in the system. The procedure by 
which this basis is constructed is the Schmidt Decomposition, which we now 
describe. 
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3.1 Schmidt Decomposition 

An arbitrary wavefunction \^) can be rewritten as a sum of two state vectors 
through a transformation of the basis vectors in the following way j^. Let \u) 
and \v) be unit vectors of the first and second particles respectively. Define M 

by 

M^{u,v\ip). (6) 

|Mp is nonnegative and bounded so it attains its maximum, for some 

and The choice of and \vi) is not unique since they can at most be 
determined only up to phase and it is possible that there are other degeneracies 
as well. 

Let \u') be any state of the first particle which is orthogonal to \ui), and let 
e be an arbitrarily small number. Then 

\\m)+e\u')f = l + 0{e') (7) 
so that up to order e^, |wi) + e\u') is a unit vector. We then find that 

{ui+eu',vi\i;) = Mi + eiu^viltP) (8) 

so 

\{ui + eu',vi\^)f = \Mi\^ + 2Re{Me{u\vi\ij))+Oie''). (9) 

But \ui) was chosen so that the scalar product in Eq. (j6|) is a maximum. There- 
fore we must have that 

\Mi\^ + 2Re(M;e{u',vim < \Mi\^, 

=^ Re(M[e{u,vi\TP)) <0. (10) 

The choice of phase of \vi) is arbitrary, however, so to ensure that Eq. (|l^) is 
satisfied we must have that 

{u',vi\^p) ^0, y u' en'^' ^{u eni\{ui\u') ^0}. (11) 

Similarly, we can show that the same restriction applies to the second particle 
so that 

{ui,v'\i;) = 0, \/ v' en"/ en2\{vi\v')=0}. (12) 

If we now define a new wavefunction 

I^A') - lV^>-Mi|ui,«i) (13) 
then \ip') also satisfies Eqs. ( pT[ ) and ( p^ ) and has the additional property that 

{uuvi\^')=0. (14) 

From this it follows that Tp' e 7i"^ (g) TC^^ and we can repeat the above process 

on ^p' . Importantly, the dimension of Ti."^ is smaller than the dimension 

of Hi (S) TL2 so this process must terminate. We finally obtain 



= ^M,\u,,v^), (15) 
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where the sum is over the smaher dimension of Tii and 7^2 and and 
are orthonormal sets. 

For the case where each subsystem has dimensionahty 2 we find that 

=Mi\ui,Vi) +M2\U2,V2). (16) 

The phases of Mi and M2 can be absorbed into and |u2,f2) so we can 

set them to be real. This decomposition can be written expHcitly as 

MX 112 ly f , ,o^ »r/2 /^cos(6'i/2)e-*'>i/2\ /cos(e'2/2)e-*''^^/2\ 
=^ ' (^cos(a/2)e-/^ sin(^;/2)e^^^/^ ) ^ \ sin(02/2)e^^^/^ ) 

+ sinfa/2)e-^/^ ^^^^^i/^)^"*^'^ « / sin(02/2)e-*/^ \ \ 

This is the Schmidt decomposition for a bipartite 2-state system. In writing this 
we have satisfied the condition that cosa > sina since otherwise |(it2, i;2|V')|^ > 
which contradicts our choice of \u\,v\). 



3.2 The Density Matrix 

If we want to calculate expectation values for one particle only and the state 
of the other particle is unknown then clearly we cannot write down the full 
wavefunction. We are forced instead to turn to the density operator, p, defined 

by 

p = \i,m= (18) 

In terms of the density operator the expectation value of any observable Q is 
given by 

(Q) - tr(pQ). (19) 
The density operator for each particle is given by 

h = tr2/5 = ^(jIpIj) P2 = trip = Y.^i\p\i) (20) 

d i 

SO that the expectation value for the ith particle can be calculated by 

{Q)^ = tr(p.Q). (21) 

For systems entangled with an (unknown) environment the density matrix rep- 
resents our ultimate state of knowledge of the system. This has important 
consequences for the interpretation of quantum mechanics. For a recent review 
of these ideas, see Paz & Zurek 0. 



4 Geometric Algebra 

Geometric algebra (GA) is essentially Clifford algebra with added geometric 
content. Since Clifford algebras are a fundamental part of the treatment of 2- 
state quantum systems (through the description of quantum spin), we expect 
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that formulating the theory in a GA framework should bring added geometric 
insight. This idea was first explored by Hestenes in a series of papers dating 
back to the sixties ^, We start by reviewing the treatment of single- 
particle systems. These are described within the GA of 3D space, denoted Gs- 
As an orthonormal basis for this we take 

1, {ak}, {lak}, 7 = ^1(72(73. (22) 

The reverse operation (which flips signs of bivectors and trivectors) is denoted 
with a tilde, and angle brackets {M)k are used to project onto the grade-/c part 
of M. For the projection onto the scalar part we simply write (M). For an 
introduction into the geometric algebra of 3D space see . 

4.1 Single-Particle Systems 

The simplest example of a 2-state system is provided by quantum spin. Spin 
states can be represented as complex 2-component vectors known as spinors. 
These can be given a more natural encoding within G3 by defining a linear 
one-to-one map between the state (as a complex vector) and a multivector con- 
structed from even-grade terms (scalars and bivectors). The simplest such map- 
ping is defined by ^ 

so that the basis elements |0) and |1) map as 

|0) ^ 1 and |1) ^ -7(72. (24) 

In this way tp sits inside the space spanned by {1, 7t7fc}, {k = 1,2,3). It follows 
that 

ipijj = (flo + aklak)iaa - atlcFk) = + al + aj + aj = ijjij = p (25) 

where p is the scalar magnitude of the state vector. The multivector -0 can then 
be written as 

^ = p^'^R. (26) 

R is then an even, normalised multivector in 3-dimensions and so is a rotor — 
a generator of rotations. 

The action of the Pauli matrices of Eq. (||) is given by 

^ cTfc Vcts = -7(7feV7(73. (27) 
It follows that multiplication by i is represented by 

= (71(72 0-3 1?/') ^ l/'7(73. (28) 
To construct observables we define the inner product of two spinors, 1}} and (p 

by 

{i^\<P) ^ (0, 4>)s = {H) - {(l>I(J3^)i. (29) 
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As we will see shortly, this definition generalises simply to multiparticle systems. 
From Eq. ( p9| ) the probability density is 

(^1^) ^ (tA, V)s - (V-V^) - (^/fT3^)/a3. (30) 

But Tplasip reverses to give minus itself, so it contains no scalar part. This 
leaves 

im ^ii^,i')s^{i^4'). (31) 

For the 1-particle case i/'V' is purely a scalar and is equal to p. In the more general 
case of n-particles we cannot assume that iprp is purely scalar and so Eq. ( |3l| ) , 
suitably normalised, provides the most general definition for the probability 
density. 

The other observable we can construct is the expectation for the spin in the 
/c-direction. This is given by 

= - p-^Iak ■ {^Ic73^)2 - p-\lc7k^4)Ia3. (32) 

Since ^ijj is a scalar [lakipi^) — 0. Also, ^lastp reverses to give minus itself and 
has even grade, so is a pure bivector (again in the multiparticle space we cannot 
make this assumption). Therefore, using Eq. ( p6| ) we can define the polarisation 
bivector by 

P = {p-'^Ia3rp)2 = {RIa3R)2 (33) 

so that 

Pk = (McrkW/iM^) - -I<Jk ■ {RIaiR)2 = -I(Tk ■ P. (34) 

In this way the spin of the particle can be thought of as a rotation of the 
Jo's plane, where the rotation is given by the wave function of the particle. 
The expectation value for the polarisation in the fc-direction is then simply the 
component of P in the fc-direction. This was Hestenes' original insight. A 
challenge is to extend these ideas to the multiparticle framework. 

As an example, consider the wavefunction employed in the Schmidt decom- 
position ( p7| ) 



/cos(0/2)e-*/2\ 
sin(0/2)e^*/2 ) ■ ^-^^^ 



In our single particle space this becomes 

iP = cos(6l/2)e-"*'^'^-^/2 _ sin(6l/2)/cr2e'^'^'"^/^ 
= e"'^-f'"^/2(cos(6l/2) - sin(6'/2)/cr2) 

= (36) 
and the polarisation vector P is simply 

= sin(6') cos(0)/cti + sin(6') sin((/))/cr2 + cos{9)Ia3 (37) 
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4.2 2-Particle Systems 



States for 2-particle systems are constructed in the multiparticle spacetime al- 
gebra (MSTA), which is built from n-particle relativistic configuration space. A 
basis for this space is provided by the vectors {7^} where the superscript labels 
the individual particle space. Vectors from different spaces are orthogonal and 
so anticommute. It follows that bivectors from different spaces commute, and 
hence the even subalgebra of the MSTA contains the tensor product of a set 
of non-relativistic algebras. This is precisely the algebra needed to construct a 
multiparticle wavefunction. A basis for a 2-particle wavefunction is provided by 
sums and products of the 1-particle basis elements {1, laj , la^}, where 

'^k^lklS, I'=l^lll2ll. (no sum). (38) 

Again, the superscript denotes the particle label, and we abbreviate I^o'\ to 
lal, etc. 

Currently our basis set gives 4 x 4 = 16 real degrees of freedom, whereas 
we should have only 8 for a 2-particle state. The solution to this problem 
is to demand a consistent meaning for the unit imaginary. In each separate 
space multiplication by the imaginary corresponds to right multiplication by 
/(T3. Since our new space has two such bivectors we require that 

^Ial=i^Ial ^^^{l-Iallal). (39) 

We therefore define the 2-particle correlator ^ 

E^\{l-I<7llal), E^^E. (40) 

-E is a projection operator and reduces the number of degrees of freedom by a 
factor of 2. The complex structure in the 2-particle algebra is now defined by 
the non-simple bivector J, where 

EIcjI - Elal - i {IcjI + lal) ^ J. (41) 

The two particle spinor jV', is now mapped to the multivector 

1^,0) ^V^i02^, (42) 

where the superscripts again denote which space the multivector inhabits. The 
result of the action of the unit imaginary becomes 

4)) ^ ^^(ti^EIal ^ i}}^(t)^EI(jl = ^?(t)'^J. (43) 

Consistency in this formulation is ensured by the results 

J2 = -E and J = JE = EJ. (44) 

The action of the individual Pauli matrices now becomes, for example, 

o-fc (g) / I^A) ^ -IcrlipJ (45) 

where / is the 2x2 identity matrix. A similar result holds for the second particle 
space. The action on the right-hand side keeps us in the space of correlated 
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products of even elements of C/3 . The quantum inner product is replaced by the 
operation 

(V-l^) ^ (V', 4>)s = 2{(j)E^) - 2{^J^)i. (46) 

The factor of E in the real part is not strictly necessary as it is always present 
in the spinors, but including it does provide a neat symmetry between the real 
and imaginary parts. The factor of 2 is included to ensure complete consistency 
with the standard quantum inner product. (In the general n-particlc case a 
factor of 2"'"'^ is required.) 

In Section we found that a general 2-particle wavefunction can be written 
in the form of Eq. (|l^) . To find the geometric algebra form of this wc first define 
the spinor 

We also need a representation of the orthogonal state to this, which is 
(-ci!(^%K'W2) - sin(e/2)e~-30/2^_,,3(^/2),,^,/.3./2 

=yj{e,c^)Ia2. (48) 

It is a straightforward exercise to confirm that this state is orthogonal to ip{9, 0), 
as required. We can now construct the MSTA version of the Schmidt decompo- 
sition. We replace Eq. ( p7| ) with 

ij (cos(a/2)i^i(0i, 0i)V;2(02, 02)e^"/' 

+ sin(a/2)V''(0i,0i)V''(02,02)/f72^fT|e^^"/')e^^S 
=pi/V'(^i,'^i)V''(^2>2)e''"/2 (cos(a/2) + sin(a/2)/a^/a^) e^'^^. (49) 
If we now define the individual rotors 

i^ = V'(0l,</'l)e^"3"/^ ^ = ^(02>2)e^"3"/^ (50) 
then the wavefunction ■0 can be written compactly as 

^P = p^/^R^S^ (cos(a/2) + sin(a/2)/cr2/cr^) e'^^E. (51) 

This gives a neat, general form for an arbitrary 2-particle state. In particular, 
all reference to the tensor product has been dropped in favour of the somewhat 
simpler geometric product. The degrees of freedom are held in an overall mag- 
nitude and phase, two separate rotors in the individual particle spaces, and a 
single entanglement angle 0. In total this gives 9 degrees of freedom, so one of 
them must be redundant. This redundancy is in the single-particle rotors. If 
we take 

R^Re^"^^, S^Se-'''''^ (52) 

then the overall wavefunction ijj is unchanged. In practice this redundancy is 
not a problem, and the form of (pll) turns out to be extremely useful. 
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The GA form of the Schmidt decomposition in is very suggestive of 
a more general pattern. To the left we have rotation operators in each of the 
individual spaces. In one sense the rotors R^S^ can be viewed as representing 
the nearest direct product (separable) state. Next comes a term describing the 
2-particle entanglement. The generalisation seems fairly clear. For a 3-particle 
system we expect to see terms describing the various 2-particle entanglements, 
followed by a term for the 3-particle entanglement. Finding precisely the optimal 
decomposition along these lines is an open problem, but the GA formalism 
has suggested an approach to the general problem of classifying multiparticle 
entanglement which has not been tried before. 

4.3 2-Particle Observables 

We can start to appreciate the utility of the form of ( pT| ) by studying the 2- 
particle observables. These go as, for example 

Wak^i IV') ^ (V-, -If^l^J). = -S/fji • (i^ JVi) (53) 

and 

(V'l a, (g, cTk 1^) ^ (V, -/a] IaliP)s = -2{I<j] laD-iijEi;). (54) 

All of the observables one can construct are therefore contained in the multivec- 
tors tpEip and tpJip. This is true in the general n-particle case, and is a major 
strength of the MSTA approach. 

To study the form of the observables we first simplify slightly and set p — 1. 
We find that (using EE = EE ^ E) 

tPEi) = R^S'^{cos{a/2) + sm{a / 2) I I al)E 
(cos(a/2) + sm{a/2)Ia^Iai)R^S^ 
= R^S^{1 + s\n{a)Iallal)ER^S^ (55) 

Substituting in the form of E from Eq. (|4^) gives 

i^E'4} = \R^S^{l - lallal + sin(a)(/cr2^/cr^ - Ia\l(jl)^R^ S'^ . (56) 
To make this result clearer we introduce the notation 

Afc = RlauR, Bk = SlakS (57) 

so that 

xPE^ = i(l - A\Bl) + \ sin(a)(^ii32 _ ^i^2)^ (gg) 
From this we see that 

m^,) = \. (59) 

This factor of one-half is absorbed by the factor of 2 in the definition of the 
quantum inner product ( p6| ) and shows that the state is correctly normalised 
to 1. The 4- vector part of the observable is more interesting, as it contains 
combinations of Ai, A2, Bi, B2, none of which are accessible to measurement 
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in the single-particle case (as they are not phase invariant). This is one place 
where differences between classical and quantum models of spin start to emerge. 

The second observable to form from the 2-particle state ip is ipJ'^Pt which is 
given by 

V'JV' = i?^S'2 (cos(a/2) + sin(Q!/2)/cr2^/(T|) J 
(cos(a/2) + sin(a/2)/cr^/cr|)^i^2 

= ^R^ S"^ (cos'^ {a/2) -sin^ {a/2)) {lal+IaDR^S^ 

= ^cos{a){Al+Bi) (60) 

This result extends the definition of the polarisation bivector to multiparticle 
systems. An immediate consequence of this definition is that the lengths of the 
bivectors are no longer fixed, but instead depend on the entanglement. 

4.4 The Density Matrix 

The density matrix for a normalised 2-particle pure state can be expanded in 
terms of products of Pauli matrices as 

p = = i(/(g) / -f- afeCTfc (g) / -f (Tfe -f Cjk (Tj «) (Tfe). (61) 

The various coefficients are found by forming, for example 

ak = {^\(7k®i\^)^-2Ial{^JiP) (62) 

It follows that all of the degrees of freedom present in the density matrix are 
contained in the multivector observables TpEip and ipJip. For mixed states we 
simply add the weighted values of these observables formed from the pure states. 
This picture is quite general and works for any number of particles. One small 
complication is that the terms in ijjj^ are anti-Hermitian, whereas the density 
matrix is Hcrmitian. One way round this is to correlate all of the pseudoscalars 
together and map all bivectors back to their dual vectors |l3| . One can often 
ignore this feature, however, and work directly with the observables ipE^ and 

An advantage of this way of encoding the density matrix is that the partial 
trace operation to form the reduced density matrix simply consists of throwing 
away any terms in the observables coming from spaces where the state is un- 
known. For example, taking the 2-particle entangled state ( |5l| ) and tracing out 
the degrees of freedom in space 2 just leaves 

p= i(l + Pfeafe), Pk = {-Iak)-{cos{a)RIa3R). (63) 

This shows that the effect of the entanglement is to reduce the expectation value 
for the polarisation from 1 to cos(q!), but leave the direction of polarisation 
unchanged. For 2-particle pure states we also see that the polarisation vector 
is the same length for both particles, so each particle is effected equally by any 
entanglement which is present. For higher particle number or mixed states the 
effects of entanglement are more complicated, though the formula 

Pfc = -2"-i(/a^).(^JV;) (64) 
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holds whenever we form the reduced density matrix for particle a from a larger, 
entangled state. 

Our simple 2-particle system exhibits one of the basic results of quantum 
theory. When a system entangles with a second, unknown system (usually 
the environment), the state of the system of interest can no longer be known 
for certain and we are forced to adopt a density matrix viewpoint. That is, 
entanglement with the environment leads to decoherencc and loss of information. 

A useful application of the preceding is to the overlap probability for the 
inner product of two states. Given two normalised states we have 

Pii^,<P)^\{mf =i<hP4')- (65) 

The degrees of freedom in the density matrices are contained in tpEip and ipJtp, 
with equivalent expressions for One can then show that the probability is 
given by the compact expression 

P(V, <!>) = {{i>Ei,){<t>E4,)) - {{i,Ji,){<f>J4,)). (66) 

This formula holds in the n-particle case as well, except for the presence of 
an additional factor of 2"~^ to give the correct normalisation. This compact 
expression is a unique feature of the MSTA approach. 

As a check on the preceding, suppose we have two separable states 

V' = R^S'^E, (j) = U^V'^E (67) 

with 

i)Ei} = ^{l-A^B^), <PE^ = ^{1-C'D^). (68) 

We find that 

P(V,,(/,) = i((l - A^B^){1 - C^D^) - {A^ + B^){C^ + D^)) 
= 1(1 + A-C B-D - A-C - B-D) 

= i(l-A-C)i(l-B-D) (69) 

which shows that the probability is the product of the separate single-particle 
probabilities. If one of the states is entangled this result no longer holds. 



4.5 Example - The Singlet State 

As a simple example of the some of the preceding ideas, consider the spin-singlet 
state 

1^) = i=(|01) - |10)) - V = -^{Icrh - IctDE. (70) 

This state is maximally entangled (a = 7r/2), and isotropic. Forming the two 
observables we find that 

i,Ei,= l{l + Iallal) (71) 

and 

iiJi) = 0. (72) 
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It follows that the reduced density matrix for either particle space is simply 
one-half of the identity matrix, and so all directions are equally likely. If we 
align our measuring apparatus along some given axis and measure the state of 
particle one then both up and down have equal probabilities of one-half. 

Suppose now that we construct a joint measurement on the singlet state. 
We can model this as the overlap probability between 'ip and the separable state 

= R^S'^E. (73) 

Denoting the spin directions by 

Rla^R ^ P, SlasS = Q, (74) 

we find that, from (|66| ) 

P{i,,cj>) = {l{l-P^Q^)\{l + Iallal)) 
= \{l-P-{Iak)Q-{Iak)) 

= 1(1- cos 0) (75) 

where is the angle between the spin bivectors P and Q. So, for example, the 
probability that both measurements result in the particles having the same spin 
{9 = 0) is zero, as expected. Similarly, if the measuring devices are aligned, the 
probability that particle one is up and particle two is down is one-half, whereas 
if there was no entanglement we should get the product of the separate single 
particle measurements (resulting in 1/4). 

It is instructive to see how all of strange quantum entanglement results for 
the singlet state are contained in the 4-vector part of the observables. This 
reveals some of the complex geometry associated with multiparticle quantum 
mechanics. And this is only for 2-particle systems! Most proposals for quantum 
computers have in mind a far greater number of entangled qubits. We hope that 
this paper has demonstrated some of the potential power of geometric algebra 
for helping to navigate through these large Hilbert spaces. 
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